Low Energy Singlets in the Heisenberg Kagome Antiferromagnet 
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The spin half Heisenberg antiferromagnet on the Kagome lattice, is mapped by Contractor Renor- 
malization to a Spin-Pseudospin Hamiltonian on the triangular superlattice. Variationally, we find 
a ground state with columnar dimer order. Dimer orientation fluctuations are described by an effec- 
tive 0(2) model at energies above an exponentially suppressed clock mass scale. Our results explain 
the large density of low energy singlets observed numerically, and the non magnetic T^ specific heat 
observed experimentally. 

PACS numbers: 75.10.Jm, 75.10.Hk, 75.30.Ds 



Frustrated quantum antiferromagncts (AFM) are im- 
portant paradigms of emerging phenomena in models of 
condensed matter. It has long been appreciated that 
due to the extensive classical ground state manifold of 
the Kagome lattice|l|, quantum fluctuations may destroy 
magnetic order and replace it with paramagnetic phases 
with novel excitations at low energy scales. 

Numerical studies of the spin half Kagome AFM have 
suggested that its ground state does not support long 
range spin order [3 S m- The low spectrum of the 
Kagome i^ consists of singlets whose number increases 
exponentially with the lattice size. 

Experimentally, there is amounting evidence of un- 
usual low energy excitations in Kagome like systems. In 
spin-| SrCrgGai20igip^, i^ a significant fraction of the 
spin moment is not frozen below the non linear suscepti- 
bility maximum at T=5K. Recently, niuon resonance ex- 
periments on a spin-i system 8], reported that below the 
susceptibility maximum of T=2GK, low frequency spin 
fluctuations were detected but without static magnetiza- 
tion down to 50mK. 

However, the specific heat of SrCrgGauOiQ has an 
unexplained large T^ coefficient which apparently does 
not arise from spin wavesQ- 

Thus there is both numerical and experimental evi- 
dence that there are seemingly non magnetic massless 
modes whose origin has not yet been understood. Weak 
bonds perturbation theory has been applied |9l[l0l|. and 
interesting results have been found for the Quantum 
Dimer Model (QDM)[1|T3 on the Kagome lattice. How- 
ever, the QDM has not yet been quantitatively derived 
from the Heisenberg model. 

It is the purpose of this Letter to elucidate the nature 
of low energy excitations of the S=l/2 Kagome AFM 
by app lyin g the Contractor Renormalization (CORE) 
method|l2l|. which is not perturbative in weakened bonds. 
CORE has been recently applied to the square lat- 
tice Hubbard model J^J, Heisenberg ladders [l4) (includ- 
ing detailed convergence tests [13), and the frustrated 
Checkerboard and Pyrochlore latticesjlq. 

For the Kagome lattice, CORE leads to an effective 
Spin-Pseudospin (S-L) model on a triangular lattice. A 
variational analysis reveals columnar dimer order in the 
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FIG. 1: CORE on the Kagome lattice (solid circles). Triangu- 
lar blocks of first (second) CORE steps are encircled by solid 
(dashed) lines. On the right: triangle four ground states are 
labelled by spin (arrows) and pseudospins (wide arrows) . The 
^■(l state of the triangle is visualized as a dimer singlet on the 
bottom rung. 



spin disordered ground state, and low energy excitations 
which can be understood as dimer orientation fluctua- 
tions. We describe these fluctuations by a p = 6 Quan- 
tum Clock Model. Its mass gap is strongly suppressed 
by quantum fluctuations by a factor estimated at about 
10^"*. Thus, the low energy singlet spectrum is in effect a 
quasi-Goldstone mode of an 0(2) order parameter. The 
number of sub gap singlets increases exponentially with 
lattice size and gives rise to a T^ specific heat as seen 
exp er iment ally. 

The Heisenberg Hamiltonian on the Kagome lattice, 
(see Fig. Pl, is 
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Henceforth we set the unit of energy to J = 1. 

Blocks. CORE involves (i) an initial choice of elemen- 
tary blocks which cover the lattice, and (ii) a truncated 
set of block eigenstates whose tensor product spans the 
reduced Hilbert space. It is useful to choose minimally 
sized blocks which respect (as much as possible) the lat- 
tice point group symmetry. Here we select upward tri- 
angles, and a truncated basis of the four degenerate spin 
half ground states, discarding the higher S — 3/2 states, 
(see Fig.©). 



The S-L representation of the four ground states are 
labelled by \s, I), where s =1 , | is the magnetization and 
I =11", -IJ- is the pseudospin in the z direction. Explicitly, 
in the Ising basis \s1S2S3), 
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The pseudospin direction in the xz plane correlates 
with the direction of the singlet bond, e.g. ■^ describes a 
singlet dimer on the bottom (— z) edge (see FigQJ. Thus, 
the L^ eigenstates have definite chiralities. 

Effective Hamiltonian. The effective interactions be- 
tween triangles is calculated by CORE[l3,[i3- We note 
that this approach is feasible when two conditions are 
met: (i) Interaction matrix elements fall off rapidly with 
range such that the truncation error at finite ranges is 
small, and (ii) the norms of the projected eigenstates are 
sufficiently large for numerical accuracy. We have com- 
puted all range 2 and range 3 interactions, and neglected 
range 4 corrections, whose expectation values were found 
to be an order of magnitude smaller. At range 3, norms of 
projected eigenstates were greater than 0.75, with most 
states above 0.9. 

The effective Hamiltonian is a Spin-Pseudospin (SL) 
Model on the triangular lattice: 
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Here L-'- — (L^, L^), and eij,e\j are unit vectors in the xz 
plane. TCss describes interactions of the Kugcl-Khomskii 
tvDe|lCl llTJ . where the pseudospin exchange anisotropy 
depends on the sites and bond directions. (See Figl^J. 
For any other other bond {ij'), Qiji is simply found by 
rotating e^ by 0,±120° according to the 0(2) rotation 
of (ij) -^ {if). 

The coupling constants and angles of 'HsL^ are tabu- 
lated in TableQ] Missing from (|3J| are terms which vanish 
in the periodic lattice by summation over nearest neigh- 
bors, and three site (ring exchange) interactions. The 
largest term KSi ■ SjLfL^L^ has a largest matrix ele- 
ment of magnitude 0.02. 

While HsL may prove to be useful for numerical stud- 
ies of the spectrum of larger lattices, its complexity some- 
what obscures its physics. It is simple however to stud 
"HsL variationally using pseudospin coherent states 
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FIG. 2: CORE range 2. Directions of anisotropy vectors e^j = 
elj = Bij for a horizontal bond. For other bond directions, 
vectors must be rotated by ±120°. The ground state singlet 
correlations of two coupled blocks are depicted by thick lines. 
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TABLE I: CORE up to range 3: Interaction parameters for 
the effective hamiltonian TisL, Eq. (|5Jl. Underlined are the 
"Dimerization Fields" (see text). 



minimized with respect to the directions li 



E^'^^[l\ = Y,{SSj)^^F[l]+Eu[l] 
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We start by evaluating the energy of the magnetically 
ordered state, where both the spins and the pseudospins 
form three sublattice (SSL) Neel order on the triangular 
lattice (and VS x ^/i order on the Kagomc) . Other can- 
didates are the dimer coverings of two triangle singlets, 
whose correlations are defined by Fig. |21 The dimer sin- 
glet states have been shown by Mila and Mambrini id] 
to span much of the low singlet spectrum in finite clus- 
ter calculations. The variational analysis highlights the 
special role of the "Dimerization Fields", Jssien Jssiei 
in (j3l, for the formation of local singlets. These terms 
cancel under summation in all uniform states defined by 
(SjSj) = const. Their significant magnitude (see Ta- 
ble ^ helps to lower the energy considerably by aligning 
li with the anisotropy vectors e^ to form singlets on cer- 



tain bonds and not others (Si • So 



This 
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is a strong argument in favor of a paramagnetic ground 
state. Consequently, Tin is crucial in selecting the true 
ground state among the multitude of dimer singlet cover- 
ings. We have found that the perfectly ordered columnar 
dimer (CD) state minimizes Hu. A local "defect" of a ro- 
tated dimer in the CD background costs a "twist" energy 
of +0.01 per site. 

In Fig. O we depict the SSL and CD states. Their 
energies per site are 



\s,l). Its energy is 



£;c£>/site = -0.229, 



EssL/site = -0.178, . (5) 




FIG. 3: Variational states on the triangular superlattice: The 
Three Sublattice (SSL) and the Columnar Dimer (CD) State, 
arrows denote pseudospins, and thick lines denote singlet 
dimers. 




FIG. 4: Clock fields <?!)(x) (thick arrows) defined by the dimer 
directions irom the sites marked by large circles. 



where the evaluation uses the known spin correlations 
of the Heisenberg AFM on the triangular lattice p^ 
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0.18. 



This result can be connected to numerical diagonal- 
izations data as follows. The number of CD states on 
Kagome clusters with the lattice group symmetry is 24: 
From a particular up-triangle site there are 6 dimer di- 
rections. There are two equivalent dimer orderings of 
the neighboring lines of dimers. Another factor of two is 
given by the down triangle configurations. Between these 
24 CD states there are exponentially vanishing overlaps 
at large lattices. 

Quantum Clock Model. A continuum Hamiltonian for 
the low energy fluctuations is derived as follows. Using 
the subset of site positions x belonging to a 2 x 2 super- 
lattice (see Fig. 0J|, every dimer configuration defines a 
unique configuration of 6-fold clock angles 0(x)'s defining 
the orientations of dimers from the selected sites. 

A ferromagnetic state of (p's represents a CD ground 
state (up to global translations, with vanishing overlap, 
of an interpolating line of dimers). The resistance to 
local twists, governed by Jieie, is described by energy 
density \psiy(j))^,Ps — 0.01. The barrier height between 
dimer orientations is estimated from Eq.ljSJ to be Iiq — 
0.05 which defines the "clock field" /ig cos(60). JiyiyLyLy 
interactions rotate the pseudospins and the dimers in the 
plane, giving rise to a kinetic energy 2X^1 1/x — 0.01. 



Thus we arrive at a long wavelength partition func- 
tion of a (2-hl)D Quantum Clock Model (QCM) which 
describes the lowest singlet sector of the Kagome spec- 
trum: 



^singlets = J ^^'"''^{Jo '^^'^''^ 



Ix'P^ + ^Ps(V0)' + he cos(6(/.) ) (6) 



The renormalization group analysis of Jose et. aL[20j 
for the classical p-state Clock Model in two dimensions 
found that at low temperatures T < Tp = Sn/p^ the 
clock field hp is a relevant interaction which locks the 
ground state into a clock minimum, with a finite gap for 
domain wall excitations. For our QCM the clock term 
is still relevant, but quantum fluctuations, which involve 
tunnelling between clock minima, drastically renormalize 
down the value of the clock field and the mass gap for 
the long wavelength excitations. 

This is shown as follows: we expand the action of © 
to lowest order in he cos(60) and integrate out the high 
momenta and energy modes 0> 



V (/)>cos(p((/)< +(/)>)) exp( --(/)>G' ^ 



= he cos{p(l)<)e-p' '''''>'>, 



hl^^ = hee 



'Cp' 



0.05-10" 



(7) 



Our rough numerical evaluation of C « 0.2 uses the 
zero point phase fluctuations of an effective spin 2 quan- 
tum xy model, describing four pseudospin half in the unit 
cell. Eq. ^ is our key result: tunnelling between p — Q 
ground states renormalizes down the clock mass gap by 
a gaussian function of p. In particular, the dispersion of 
6 fluctuations 



•'=' = j{hrY + '-\ce 
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resembles Goldstone modes of an 0(2) xy-model at fre- 
quencies and temperatures larger than Q. In this 
regime, the singlets' contribution to the specific heat is 
quite large. Using thermodynamics of free bosons, we 
obtain 



C\ 



—T^, S{E) 

Ps 



N 



E 



7V0.01 



2/3 



(9) 



where N is the number of effective sites. Thus, the small- 
ness of /ig*^" provides the long sought after explanation 
of the unusual exponential proliferation of singlets at 
sub gap energies j5|. In addition, the singlets pseudo- 
goldstone mode can now explain the experimentally re- 
ported T^ term in the specific heatQ. 

We note that the CD state has no long range spin or- 
der. The spin gap can be estimated from the variational 
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0.113 0.08 -0.005 0.026 0.182 -0.039 330° 280° 
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TABLE II: Second CORE iteration: Interaction parameters 
for TisL evaluated up to range 2. 



lattice size 3 9 27 81 



spin gap 1 0.5 0.1 0.06 



TABLE III: Spin gap as a function of effective Kagome lattice 
size using the Heisenberg model, first and second iterations of 
CORE. 



energy difference between CD and SSL states to be of the 
order of hg = 0.05. Another estimate can be obtained by 
iterating CORE on Hsl- 

At the second CORE step, the triangular lattice is 
covered by triangles (which are blocks of 9 Kagome lat- 
tice sites) . These form a triangular superlattice with di- 
rected bonds (see Fig^. For both the Heisenberg and 
Ti-SL, each block has four degenerate S = 1/2 ground 
states which can again be represented by a spin and a 
pseudospin. The second CORE step thus maps the S-L 
Hamiltonian Q onto itself with new interaction param- 
eters and anisotropy vectors, as listed in Table ITU 

What can we learn from step 2? 

In contrast to the first CORE step (see Tabled, the 
vectors e^j • eji > 0, i.e. ferromagnetic. Jssieie has de- 
creased while Jss did not. Thus the Hamiltonian prefers 
local singlet correlations with aligned pseudospins, which 



is consistent with columnar order. 

Spin gap. By iterating CORE we can compute the 
splitting between the S — ^ ground state and the lowest 



S — ^ excitation on triangular clusters. In Table IIIII 
the spin gap is computed on up to 81 original Kagome 
lattice sites. At third CORE step, many wave function 
overlaps vanish. This is expected due to the onset of 
long range CD order, since states with high pseudospin 
parentage have lower energy. Stopping at 81 sites, we 
can only observe that while the spin gap is larger than 
the dimer fluctuations bandwidth, there is no conclusive 
sign of saturation to a finite thermodynamic limit. 

In summary, we have computed the effective Hamilto- 
nians of the S=l/2 Heisenberg antiferromagnet on the 
Kagome lattices using CORE up to range 3. Variation- 
ally, we conclude that the Kagome lattice has long range 
columnar dimer order and a very low gap for singlet ex- 
citations in the thermodynamic limit. CD order might 
induce experimentally detectable lattice distortions. Fur- 
ther details can be found else where l2ll. 



Note added: In a recent preprint, Capponi, Laeuchli, 
and Mambrini|23| have independently computed TCsl of 
Eq. (PI by CORE, and found excellent numerical agree- 
ment with large cluster diagonalizations. 
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